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Abstract. The Kac-Ward formula allows to compute the Ising parti- 
tion function on a planar graph G with straight edges from the deter- 
minant of a matrix of size 2N, where TV denotes the number of edges 
of G. In this paper, we extend this formula to any finite graph: the 
partition function can be written as an alternating sum of the determi- 
nants of 2 2fl matrices of size 2N, where g is the genus of an orientable 
surface in which G embeds. We give two proofs of this generalized for- 
mula. The first one is purely combinatorial, while the second relies on 
the Fisher-Kasteleyn reduction of the Ising model to the dimer model, 
and on geometric techniques. As a consequence of this second proof, we 
also obtain the following fact: the Kac-Ward and the Fisher-Kasteleyn 
methods to solve the Ising model are one and the same. 



1. Introduction 

The "algebraic method" discovered in the 40's by Onsager and Kauf- 
man [19[ I16j to solve the Ising model on a square lattice G was widely 
considered to be extremely difficult. This motivated Kac and Ward [12| 
(and later, Potts and Ward |20j) to try to find a more direct "combina- 
torial method", building on van der Waerden's reformulation of the Ising 
model |24| . They defined a matrix M(G), with rows and columns indexed 
by oriented edges of G, whose determinant coincides with the square of 
Z I {G), the Ising partition function on G. Unfortunately, many arguments 
in \12\ [20] are of heuristic nature, and some key topological statement later 
turned out not to hold [21]. 

Since the technicalities raised by a rigorous proof of the Kac-Ward for- 
mula seemed formidable, the focus shifted to finding combinatorial methods 
not involving directly the Kac-Ward matrix M(G). This was achieved in- 
dependently and almost simultaneously by Sherman [21], Hurst-Green |10j . 
Kasteleyn |13] and Fisher [6 J . Sherman built on (unpublished) ideas of Feyn- 
man to express Z T (G) as a formal infinite product of polynomials, for any 
planar graph G with vertices of degree 2 or 4. The other three authors 
all essentially obtained the following result: they related Z T (G) with the 
dimer model partition function Z d (Tq) on an associated graph Fq, and then 
found a skew-symmetric adjacency matrix A(Tq) for Tq whose Pfaffian gives 
Z d (Tg)- This "Pfaffian method" was later extended by Kasteleyn [14"1 [T5] 
to any planar graph. 



2000 Mathematics Subject Classification. 82B20, 57M15, 05C90. 

Key words and phrases. Ising model, Kac-Ward matrices, spin structure, dimer model, 
Kasteleyn matrices. 

1 



2 



David Cimasoni 



To the best of our knowledge, the first direct combinatorial proof of 
the Kac-Ward formula Z 1 (G) 2 = det(M(G)) for any planar graph G with 
straight edges was only obtained in 1999 by Dolbilin et al. [5]. 

In this paper, we present a generalized Kac-Ward formula, valid for any 
finite graph G (possibly non planar, disconnected, with loops and multiple 
edges). More precisely, let £ be an oriented closed surface of genus g in 
which G embeds. For each spin structure A £ S(£), we define a generalized 
Kac-Ward matrix M X (G) whose rows and columns are indexed by oriented 
edges of G (see Definition [1]) . Our formula reads 

Zl ^ = h £ (-l) Arf ( A ) det(M A (G)) 1 / 2 , 
Aes(s) 

where Arf(A) 6 Z2 is the Arf invariant of the spin structure A (Theorem l2.ip . 
If G is a planar graph with straight edges, this formula is precisely the clas- 
sical Kac-Ward formula Z\G) 2 = det(M(G)). Using Bass's theorem on the 
Ihara-Selberg zeta function for graphs [I], we immediately obtain a gener- 
alized Feynman-Sherman type formula expressing Z 1 (G) as an alternating 
sum of 2 2g formal infinite products (Corollary 12. 2p . This slightly extends a 
result of Loebl [T7] (from the case with all vertices of degree 2 or 4 to the 
general case), and expresses it in what we believe to be the right language 
(spin structures instead of Sherman rotation numbers). 

Our generalized Kac-Ward formula is new and original, but its mere ex- 
istence is by no mean a surprise. First, it is known since Dolbilin et al. 
Tesler [22], Gallucio-Loebl [9] and Cimasoni- Reshetikhin [3] that the dimer 
partition function Z D (T) on any finite graph T can be expressed as the 
same alternating sum of 2 2g Pfaffians of skew-symmetric adjacency matrices 
A X (T); by the correspondence G h-> Tg, such a formula holds for Z 1 (G) as 
well. Also, as mentioned above, an analogous Feynman-Sherman type for- 
mula was already obtained by Loebl [17J. Finally, Masbaum informed the 
author that he independently obtained Corollary 12.21 the improved version 
of this Feynman-Sherman formula where spin structures replace Sherman 
rotation numbers. He was therefore aware (via Bass's Theorem) of our for- 
mula. 

However, we would like to emphasize the fact that we give two completely 
independent proofs of our main result. The first one is of combinatorial na- 
ture, following the line of Dolbilin et al. [5J. (These authors developed most 
of the necessary tools to obtain the general formula, but without the under- 
standing of the importance of spin structures, one cannot define the right 
generalized Kac-Ward matrices.) The second proof is of geometric nature. 
We first use the Fisher correspondence G 1— > Tq together with the geometric 
treatment of the dimer model [3] to obtain the formula for Z 1 '(G) in terms of 
the Pfaffians of 2 2 9 Kasteleyn matrices A X (T G ). (Th is was done in an anal- 
ogous way in [3].) Then, we explicitly transform these Kasteleyn matrices 
into the corresponding generalized Kac-Ward matrices M X (G), using moves 
that do not change the determinant. Therefore, not only do we general- 
ize the Kac-Ward approach and clarify the Feynman-Sherman one, but we 
also unify these two approaches with the Fisher-Kasteleyn Pfaffian method 
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of |1CH fl"3j El H] . The fact that the Kac-Ward and Pfaffian methods are one 
and the same is by no mean obvious, and they are very often considered as 
distinct (see e.g [23 Q7]). 

The paper is organized as follows. In Section [21 we introduce the general- 
ized Kac-Ward matrices (Definition [1]) , state our main result (Theorem 12. ip 
and its consequence (Corollary 12. 2p . In Section [3l we give the combinato- 
rial proof, following the line of [5], and assuming many of their results. (In 
this sense, the combinatorial proof given here is not entirely self-contained.) 
Finally, Section 0] deals with the geometric proof. We begin by giving an 
improved version of the Fisher construction G i— > Tq (Subsection 14. ip . We 
then recall the geometric approach to the dimer model (Subsection I4.2p and 
finally show how to identify the Kac-Ward and Fisher-Kasteleyn matrices 
(Subsection I4.3P . 

Acknowledgments. The author would like to express his thanks to Martin 
Loebl and Gregor Masbaum for stimulating discussions. 

2. Statement of the formula 

Let G be a finite graph, possibly disconnected, possibly with loops and 
multiple edges. Let us associate to each edge e G E(G) a formal variable x e . 
As discovered by van der Waerden [23], the partition function for the Ising 
model on G with zero magnetic field can be put in the form 

z\g) =j2 n 

G'CG eeE(G') 

where the sum is over all subgraphs G' of G such that each vertex of G is 
met by an even number of edges of G' . In other words, the sum is over the 
set 

Z x (G; Z 2 ) = {£ = £ e £ e e € C 1 (G; Z 2 ) | e>£ = € C (G; Z 2 ) } 

of 1-cycles modulo 2 in G. Therefore, if we denote by |£| the set of edges of 
G such that £ e is equal to 1, the Ising partition function on G can be written 

z\G) = where = n x - 

The aim of this section is to state a closed formula for this partition func- 
tion in terms of the determinant of 2 2rj "generalized Kac-Ward matrices", 
where g denotes the genus of a surface where the graph G embeds. These 
matrices will be defined using spin structures. We shall therefore start by 
recalling the main features of these geometric objects. 

Let £ be a closed oriented smooth surface endowed with a Riemann- 
ian metric, and let FY, — > S denote its orthogonal frame bundle. A spin 
structure on E is a cohomology class A £ i^ 1 (i ? S;Z 2 ) whose restriction 
to each fiber F x of this principal 5 1 -bundle is the non-trivial element of 
i7 1 (i ? 1 .;Z 2 ) = Z 2 . In other words, a spin structure assigns to each framed 
curve in E an integer modulo 2, in such a way that the following two con- 
ditions are satisfied: the boundary of any disc with constant framing has 
value 0, and the boundary of any disc with framing induced by an outgoing 
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vector field has value 1. In particular, a spin structure allows to compute a 
modulo 2 winding number for any closed curve 7 in £: simply consider the 
value assigned to 7 with framing induced by the tangent vector field along 
7. We shall denote this number by A(t*) E Z2. 

One easily checks that the set §(£) of spin structures on S is an affine 
Z2)-space. In particular, there are exactly 2 2g distinct spin structures 
on a closed oriented surface £ of genus g. A convenient way to give oneself 
a spin structure on £ is via a vector field with isolated zeroes of even index: 
the associated spin structure will be obtained by computing the modulo 2 
winding number of this vector field along a closed framed curve, with respect 
to the framing of this curve. The evenness of the indices ensures that both 
conditions described above are satisfied. 

Recall finally that spin structures on a surface X come in two classes, 
determined by their "Arf invariant": the even spin structures, with Arf 
invariant 0, and the odd ones, with Arf invariant 1 (see the definition in 
Subsection 13. 1 j) . These two classes coincide with the orbits of the action of 
the diffeomorphism group of £ on S(£). 

We are now ready to define the generalized Kac-Ward matrices. 

Let G be a finite graph, and let E denote the set of oriented edges of 
G. For any e G E, we shall denote by s(e) its starting point, by /(e) its 
endpoint, and by — e the same edge with the opposite orientation. By abuse 
of notation, we shall denote by x e = X- e the formal variable associated to 
the unoriented edge corresponding to e and — e. Now, embed G in a closed 
oriented surface £, and fix a spin structure A G §(£) via a vector field on £ 
with isolated zeroes of even index in S \ G. Finally, fix an arbitrary point 
in each edge of G C £. 

Definition 1. The generalized Kac-Ward matrix associated to G and A is 
the x \E\ matrix I — T A (G), where I is the identity and T^(G) denotes 
the transition matrix defined by 



Here, a\(e, e') E R is the rotation number of the tangent vector field along 
e followed by e' with respect to the vector field A, from the fixed point in 
the edge e to the fixed point in the edge e'. 

Obviously, the matrix / — T X (G) will depend on an ordering of the set E, 
on the choice of the vector field representing the spin structure A, and on 
the choice of the fixed points in the edges of G. However, its determinant 
will not. This is the first part of our main result. 

Theorem 2.1. Let G be a finite graph embedded in an orientable closed 
surface £ of genus g. For any spin structure A € §(£), the determinant of 
the matrix I—T X (G) is the square of a polynomial in the variables {x e } e ^E(G) 
which only depends on A and on G. If det(I — T A (G)) 1 / 2 denotes the square 
root with constant coefficient equal to +1, then the Ising partition function 




exp (|a\ (e, e')) x e if /(e) = s(e') but e' 7^ — e; 
otherwise. 
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on G is equal to 



Z ^ G ) = h E (-l) Arm det(I-T\G))^, 

AGS(S) 

where Arf(A) G Z2 is the Arf invariant of the spin structure A. 

A remarkable feature of this formula is that the left-hand side only de- 
pends on the abstract graph G, while the right-hand side a priori depends 
on the way this graph is embedded in a surface. 

Example. Let G be a planar graph with straight edges. Obviously, there is 
a single spin structure A in the plane, whose most natural representative is 
a constant vector field. The associated rotation number a\(e,e') is simply 
given by the oriented angle <(e, e') between the oriented edges e and e'. 
Therefore, we have the formula 

Z I (G) = det(M(G)) 1 / 2 , 

where M(G) is the \E\ x \E\ matrix with coefficients 

' 1 if e = 

M(G) efi , = \ - exp (|<(e, e')) x e if f(e 

otherwise. 

This is precisely the Kac-Ward formula [12J. 



e'; 



s(e') but e' 7^ — e; 



Example. Let G be the "figure-eight graph" given by a single vertex and 
two loops e±,e2 at this vertex with associated variables x\,X2- Obviously, 
the Ising partition function is given by Z I {G) = 1 + x\ + X2 + x\X2- We will 
now check that the formula in Theorem 12.11 gives the right result. 



e-2 



ei 



Let us first consider G embedded in the plane as 
illustrated opposite, with the marked fixed points in its 
edges. Let A be the unique spin structure on the plane, 
given by a constant vector field. The corresponding rotation numbers are 
given by 




a\{e2,e 2 ) = -a A (ei,ei) = 2tt, 

If the set E is ordered by {ei,e2, 
Ward matrix is given by 

/l + Xl 

I-T X {G) = 



-X2 





a\(ei,e 2 ) = a\(e 2 ,ei) = 0, etc... 
-ei,— ea}, then the corresponding Kac- 



-ixi \ 




\ ~iX2 



-x\ 





1 +X 2 


iX2 


ix\ 


l + Xl 





-X2 



-xi 

1 + X2J 



Its determinant being equal to (l + 2;i) 2 (l-|-2;2) 2 , w e get the claimed equality 
det(J-T A (G))V2 = Z I (G). 
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Let us now consider G embedded in the torus T 2 as 
illustrated opposite. Let A^i denote the vector field on 
T 2 induced by a constant vector field on the plane, and 
let \\-\ (resp. A-1,1, be obtained from Ai,i via 

a full twist along the closed curve 71 (resp. 72, 71 and 
72). Using the notation 00 := exp(i7r/4), we obtain 









72 

> — 1 






ei 

e 2 > 


1 — »- 

i 

7i 





det(I - T x ^2 (G)) 



1 — E\Xi —£2UX\ — UJX\ 

—E1LUX2 1 — £2^2 —OJX2 

—£162^1 1 — £\X\ —£\ojxi 

-£l£ 2 U;X2 -£ 2 WX 2 1 - £2^2 

2 



= (1 - £1X1 - £ 2 X 2 - £l£2XiX 2 

Using the definition of the Arf invariant (see Subsection 13. ip , one gets 



(_l)Arf(A £1 , £2 ) = 1(1 - £l 



£2 - £l£2j 



-1 if £1 = £2 = 1; 
1 otherwise. 



Hence, the alternating sum in Theorem 12.11 is again equal to Z 1 (G) as 
claimed. 



We now turn to a consequence of Theorem 12.11 

Let 7 denote an oriented closed path on a graph G. We will say that 7 is 
reduced if it never backtracks, that is, if no oriented edge e is immediately 
followed by the oriented edge — e. The oriented closed path 7 will be called 
prime if, when viewed as a cyclic word, it cannot be expressed as the product 
5 r of a given closed path 5 for any r > 2. Note that the oriented closed path 
7 is reduced (resp. prime) if and only if —7 is. It therefore makes sense 
to talk about prime reduced unoriented closed paths. Finally, recall that a 
spin structure A € §(£) allows to compute a modulo 2 winding number for 
an oriented closed curve 7 in S, which we denote by A (7) G 7Li- Note also 
that 7 and —7 have the same winding numbers. 

Applying Bass's Theorem [HE] to our generalized Kac-Ward matrices, we 
immediately obtain the following corollary, which extends results of Sher- 
man [21] and Loebl [17] . 

Corollary 2.2. Let G be a finite graph embedded in an orientable closed 
surface £ of genus g. Then the Ising partition function on G is equal to 

Zl ^ = h E (-i) Arf(A) []( 1 -(- 1 ) A( %)). 

Ags(s) 7 g:r 

where 31 denotes the (infinite) set of prime reduced unoriented closed paths 
on the graph G, and x('j) stands for Y\ e€ ^x e . □ 

Gregor Masbaum informed the author that he independently obtained 
this corollary working directly with infinite products as in [21[ 117] . Via 
Bass's Theorem, he therefore has yet another proof of Theorem 12.11 
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3. The combinatorial proof 

The aim of this section is to give a first proof of Theorem 12.11 of com- 
binatorial nature. More precisely, we first use Johnson's theorem on spin 
structures and some standard tricks to reduce our formula to a combina- 
torial statement (Proposition 13. lj ). We then build on some results of [5] to 
prove this combinatorial statement. 

3.1. Computing the determinants using Johnson's theorem. We 

shall start by giving a proof of the first part of Theorem \2.1\ that is, that 
the determinant of I — T X (G) only depends on A and G. As we will see, it 
follows quite easily from a result of Johnson [TT], that we now recall. 

As above, let £ be a closed oriented surface of genus g. A quadratic form 
on i/i(£;Z2) is a map q: HiiT,;^) —> Z2 such that q{x + y) = q{x) + q{y) + 
x ■ y for all i,t/£ Hi(Y>; Z2), where • denotes the intersection form. The set 
Q(S) of such quadratic forms is clearly an affine H 1 ^; Z^-space. Given a 
vector field A representing a spin structure, let q\ : HiCE;^) —> <^2 denote 
the map defined as follows. Represent a class a G i2i(£; Z2) by a collection 
of disjoint oriented regular simple closed curves C\, . . . ,C m in £ avoiding 
the zeroes of the vector field A, and set 

m 

q x (a) = £(A(Ci) + 1) (mod 2) 

i=l 

where A(Cj) denotes the winding number of the tangential vector field along 
Ci with respect to the vector field A. Johnson's theorem asserts that this 
assignment A 1— > q\ defines an i? 1 (S; Z2)-equivariant bijection between the 
set S(S) of spin structures on S and Q(S). The Arf invariant of a spin 
structure A G S(S) is then defined as the Arf invariant of the corresponding 
quadratic form q\, that is, the modulo 2 integer Arf (A) G I12 given by the 
equality 

( _ 1) Arf(A) = ^ £ ( _ 1)flx( a) i 
aeHi(E;Z 2 ) 

We shall need a single property of this invariant, namely that it satisfies the 
equality 

(1) J_ £ (_l)Arf(A) ( _ 1)9 x(a) = 1 

A£S(E) 

for any a G i?i(S;Z2). We refer to |18} Lemma 2.10] for the easy proof. 
Let us start the computation of the determinant 

det(J - T X (G)) =: detM A = JJ M^ {e) . 

aeS(E) eeE 

Each permutation a G S(E) decomposes into disjoint cycles, inducing a 
partition E = \_\- Ej (a) of the set of oriented edges into orbits of length 
£j{o~)- The corresponding contribution in the determinant is 

n M Ue) = n(-i)' j(CT)+i n 

edE j eG£j(cr) 
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Note that M A e is equal to 1 if e is not a loop. Therefore, the non-loop 
oriented edges that are fixed by a do not contribute to this product, and 
the corresponding orbits can be removed. By definition of M x , a permuta- 
tion a £ S(E) will have a non-zero contribution if and only if each of the 
remaining orbits Ej(a) forms an oriented closed path jj(cr) in G, such that 
the corresponding set of closed paths 7(0") = {7j(o")}j covers all oriented 
loops of G, and such that 7(0") is a so-called admissible path, that is: 

(i) it passes through each edge of G at most twice, and if so, in opposite 
directions; 

(ii) it never backtracks, that is, no 7j(cr) contains two successive oriented 
edges of the form e and — e. 

The corresponding non-zero contributions in det M A will be of the form 
IJ(-l)4M+i Yl (^. ((7)il -exp(la A ( e ,a(e)))x e ), 

j eeEj((j) 

where <^-(o-),i is the Kronecker symbol. This is equal to 

n(^»:i~ II exp (§a A (e,o-(e))) x e ). 

Let T(G) denote the set of admissible closed paths 7 = {7j}j in G, and 
T(G) the set of admissible closed paths that cover all oriented loops of G. 
We have: 

detM A = E II (V " II ex P (^x(e,a(e)))c 

7G?(G) 3 

= E n( <5< i. 1_exp (5 aA (T;)) a; (7j)), 
7 er(G) 3 

where x(jj) = Tieejj x e an d a x(lj) 1S the rotation number of the tangent 
vector field along jj with respect to the vector field A. Using the notation 
= - ex P(§«A(7j)), we get: 

detM A = £ l[(5 ejA +z( 7j )x( 7j )) 

7 ef(G) J 

= e n a+^M7,)) n ^-mt,) 

76? (G) 01^=1} 

= e ( e n^)^^)) n <iMn) 

7G?(G) Sc{fc|/ fc =l}jeS 

= E H^Mij) 

ier(G) j 

= : E ^(tMtO- 

7er(G) 
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We now claim that for any given admissible closed path 7 = {71, . . . , 7 m } 
in G, we have the equality 

m 

<1) ■■= n- eX P(§«A(7,)) = (-l)«A(7)+*(7) j 

i=i 

where q\ is the quadratic form associated to the spin structure A, and 4(7) 
denotes the number of double points of the path 7 perturbed in general 
position. (Note that this number is well-defined, as 7 never backtracks.) 
Both sides of this equality are unchanged when the path is perturbed, so it 
is enough to prove the corresponding equality for a family a = {a±, . . . , a m } 
of oriented closed curves in general position in S. 

This can be done by induction on the number 

t(ot) of double points of a. If t(a) is zero, a 

consists of a disjoint union of oriented simple Q 
closed curves. In this case, the equality above is 
precisely Johnson's definition of the quadratic form q\. On the other hand, 
smoothing out one double point of a as illustrated above transforms it into a 
new family a' of m' = m± 1 oriented closed curves with q\(a!) = q\(a) and 
t(a') = t(a) — 1. One easily checks that the left hand side of the equation 
displayed above will also be multiplied by —1, thus proving the claim. 

Summing up, we have proved the following: For any spin structure A 6 
S(E), the determinant of the corresponding generalized Kac-Ward matrix is 
equal to 

det(7 — T X (G)) = (-1) 9a(7)+ * (7) x(7). 
7 er(G) 

In particular, this determinant only depends on A € §(S) and on G. 

Given X,Y two subgraphs of G with no edge in common, let T(X,Y) 
denote the set of admissible paths that pass once on the edges of X, twice on 
the edges of Y, and that do not pass on the edges of G \ {X U Y) . Obviously, 
this set is empty unless X = |£| for some £ G Z^G;^)- Therefore: 

det(I — T X (G)) = (-1) 9a(7)+ ' (7) ^(7) 
7 er(G) 

= E E E (-ir (7)+t(7 MieiMn 2 

£eZ!(G;Z 2 ) y 7 er(|£|,y) 

= E (-i) 9A(?) E( E (-i)' (7) >(iei)*(y) 2 , 

£eZi(G;Z 2 ) y 7 er(K|,y) 

where the second sum is over all subgraphs Y of G with no edge in common 
with |£|. On the other hand, let Z\(G) denote the following twisted partition 
function 

Z X (G):= £ (-l)»<€)x(|f|). 

?eZi(G;Z 2 ) 
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Its square is equal to 

6,6 

= £ (_l)«x(€>£( £ (-i)^^)x(|e|)x(y) 2 , 

where the second sum is over all subgraphs Y of G with no edge in common 
with |£|. Therefore, the combinatorial proof of Theorem 12.11 boils down to 
the following: 

Proposition 3.1. For any £ £ Z^G;^) and any subgraph Y of G with no 
edge in common with |£|, we have the equality 

E (-!) 4(7) = E (-i)* 1 *- 

7er(|f|.y) «i+« 2 =« 

l«iln|« a |=y 

Indeed, this proposition implies the equality between Z\{G) and the 
square root det(J — T A (G)) 1//2 whose constant coefficient is equal to +1. 
It follows: 

z\g) = e 
= E £ 

= E (i E (-i) Arf(A) (-i)- w ) E 

oeHi(S;Z 2 ) AeS(E) 

= h E c- 1 )^ E (-i^^ E 

AeS(E) aGi/i(E;Z 2 ) 5;K]=" 

1 £ ( _l)Arf(A) ZA(G) 



29 

AeS(E) 



1 £ (-l) Arf Wdet(7-T A (G)) 1 /2. 



29 

AeS(E) 



3.2. Proof of Proposition 13.11 We shall now give a proof of Proposi- 
tion 13.11 building on and assuming results of [5] . We would like to insist on 
the fact that, although our proof will appear to be short, this combinato- 
rial statement is by no mean trivial: alone, the complete proof of a special 
case of Lemma 13.31 below takes 12 pages in [5], constituting the core of that 
paper. 

We shall actually evaluate independently both sides of the equality in 
Proposition 13.11 The right-hand side does not require too much effort. 
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Lemma 3.2. Fix a 1-cycle £ G Z\{G\ r L<i), a subgraph Y of G with no 
in common with |£|, and set 

z& Y) = {(6,6) g z x {G- z 2 ) 2 |6 + 6 = 6 161 n |6I = ^} • 

Then, the integer 

E 

(£x,6)ez(e,y) 

is equal to zero unless the following two conditions are satisfied, in which 
case it is equal to the cardinality of Hi(\£\; Z 2 ): 

(i) there exists c G Ci(|^| ; ^2) suc/i thatdc = dcy , where cy G Ci(G;Z 2 ) 
is i/ie 1- chain with \cy\ = Y; 

(ii) for any j3 G i?i(|£|;Z 2 ) t/ie intersection number (3 ■ £ is et/uaZ to zero. 

Proof. First note that any element (6,6) i* 1 Z(6Y) provides a 1-chain 
c = 6 + c y S Ci(|£|;Z 2 ) such that 9c = <9cy. Therefore, if there is no 
such 1-chain, then the set Z(6 Y) is empty and the integer considered equal 
to zero. If there is such a chain on the other hand, then the assignment 
(6,6) ^ 6 + °Y defines a bijection from Z(6 Y) onto the set of such 1- 
chains, with inverse c 1— > (c + cy , c + cy + £). This set being an affine space 
over 

Ker(<9: Ci(|£|; Z 2 ) -+ C (|£|; Z 2 )) = Zi(|£|; Z 2 ) = fftflfl; Z 2 ), 
we get the equality 

^ _ ^ (_l)(c+Cy+/3)-(c+ C y+/3+0 

(_i)(°+<*+0K 

/3eHl(|5|;Z 2 ) 

= E (-d*. 

as c + cy is an element of Z 2 ). This integer vanishes unless the linear 

form (3 1 — y is identically zero, in which case it is equal to |-£/i(|£|; Z 2 )|. □ 

Proposition 13.11 (and therefore, Theorem I2.ip now follows from one more 
lemma. Once again, we shall not give here a self-contained proof of this final 
result, but assume statements of [5]. 

Lemma 3.3. Fix a 1-cycle £ G Zi(G;Z 2 ), a subgraph Y of G with no edge 
in common with and let r(|£|, Y) denote the set of admissible paths in 
G that pass once on |£|, twice on Y, and not on G \ (|£| U Y). Then, the 
integer 

s(\h\,y)= e (- 1 )* (7) 

is equal to zero unless the conditions (i) and (ii) of Lemma 1 3. 2\ are satisfied, 
in which case it is equal to the cardinality o/iTi(|£|;Z 2 ). 



12 



David Cimasoni 



Proof. Following the terminology of [5j Sections 4-5], let us call an orienta- 
tion of the edges of |£| "regular" if around each vertex v of |£|, the oriented 
edges belonging to |£| that enter v strictly alternate with the edges leaving 
it. Observe that |£| admits a regular orientation if and only if the condition 
(ii) of Lemma 13.21 is satisfied. Therefore, if condition (ii) is not satisfied, 
then |£| does not admit any regular orientation, and [5j Corollary 2] implies 
that S(\£\, Y) is equal to zero. On the other hand, if condition (ii) is satis- 
fied, then |£| admits a regular orientation. Since E is orientable, one easily 
checks that this regular orientation can be chosen so that it extends to a 
"principal orientation" on |£| L)Y. (This is an orientation of the edges of |£| 
and of the doubled edges of Y so that, around each vertex v of |£| U Y, the 
oriented edges belonging to |£| U Y - both simple and doubled - that enter 
v strictly alternate with the edges leaving it.) Then, the proof of [5J Theo- 
rem 4] extends verbatim, leading to: S^l^y) vanishes unless condition (i) 
holds, in which case S(\£\,Y) = |iJi(|£|;Z 2 )|. □ 

4. The geometric proof 

We now give a second proof of our formula. This proof is in three steps. 
We first recall the definition of the dimer model, and state (an improved ver- 
sion of) the Fisher correspondence [7] between the Ising and dimer models 
(Subsection l4.ip . We then recall Kasteleyn's theory, as well as the geometric 
treatment of the dimer model initiated by Reshetikhin and the author in [3] 
(Subsection I4.2p . Finally, we show in Subsection 14.31 that the 2 2g general- 
ized Kac-Ward matrices can be naturally identified with the 2 2g Kasteleyn 
matrices for the corresponding dimer model. As a direct consequence, we 
obtain not only a purely geometric proof of Theorem 12.11 but also the fact 
that both methods - the Kac-Ward on one hand, the Fisher-Kasteleyn on 
the other - should really be considered as one single method to solve the 
Ising model. 

4.1. Reduction of the Ising model to the dimer model. Let T be a 

finite graph. A dimer covering (or perfect matching) on T is a choice of 
edges of T, called dimer s, such that each vertex of V is adjacent to exactly 
one of these edges. We shall denote by M(r) the set of perfect matchings 
on r. Let us associate to each edge e £ E(G) a formal variable x e . The 
partition function for the dimer model on V is by definition 

Z D (T) = X ( M )> where X ( M ) = II Xe ' 

MeM(r) eeM 

As discovered by Kasteleyn [13] and Fisher [61 [7], the Ising model can 
be interpreted as a special type of dimer model. More precisely, let G 
be a finite graph with variables x e associated to its edges. For each vertex 
v G V(G), choose a linear ordering of the adjacent edges. Then, consider the 
associated graph Tq obtained from G by blowing up each vertex v S V(G) 
into a cluster, as illustrated in Figure [TJ 

In this way, each vertex v £ V(G) of degree n gives rise to 2n vertices 
in Tq. Furthermore, there is a natural inclusion of E{G) into E(Tq). To 
each edge of Tq coming from an edge e of G, assign the same variable x e . 
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Figure 1. The neighborhood of a vertex v G V{G) of degree 
n, and the associated cluster in T^. 

To each other (new) edge of Tq, assign the weight 1. Obviously, the graph 
Tq depends on the choice of the linear orderings of the edges around each 
vertex. However, for any such choice, we have the following result. 

Proposition 4.1. The restriction map 

<p: C 1 (r G ;Z 2 )^C 1 (G;Z 2 ), £ e e ^ Yl ^ e 

e€E(F G ) e£E(G) 

induces a bijection M.(Fg) — > Z^G;^) so that x(\ip(M)\) = x(M) for any 
perfect matching M G JAi^o)- In particular, Z 1 (G) = Z d (Tg)- 

Proof. First note that for any M G "M(Tg), (p(M) is a 1-cycle modulo 2. 
Indeed, let v G V(G) be a fixed vertex and let n denote its degree. In and 
around the cluster coming from v, M matches exactly 2n + 5 m (v) vertices, 
where 5m{v) denotes the degree of v in the subgraph |y(M)| C G. As M 
is a matching, it matches an even number of vertices, so 5m (v) is even and 
9?(M) is a 1-cycle modulo 2. 

Then, one easily checks by induction on the degree 
of v that any £ G C\{G\ r L'2) with <9£ even at v (resp. 
odd at v) uniquely extends to a perfect matching on 
the cluster corresponding to v (resp. on the graph 
illustrated opposite). In particular, any £ G Z\{G;1j-i) 
extends uniquely to a perfect matching G M(r<3). 
By construction, (p(M^) = £ for any £ G Zi(G;7*2), so ip induces a bijection 
as claimed. The equality x(\ip(M)\) = x(M) follows from the definitions, 
and the proposition is proved. □ 




FIGURE 2. The preferred perfect matching Mq on Tq- 




14 



David Cimasoni 



In particular, there is a preferred perfect matching Mq on Tq, namely the 
unique one such that <p(Mo) = 0. It is illustrated in Figure [2j 

Note also that if the first graph G is embedded in a surface X, and if the 
chosen linear orderings around the vertices are compatible with the cyclic 
orderings induced by the embedding GcS, then the graph Tq also embeds 
in the same surface S. Furthermore, the bijection (p can easily be seen to 
satisfy the following additional property: for any M 6 M.(Tq), the 1-cycles 
M+Mq g Zi(Tg; Z2) and f(M) £ Z\(G\ Z2) define the same homology class 
in H\(^lj\!L<i) . Therefore, for any a € i?i(S;Z2), we have equality between 
the corresponding partial partition functions for the Ising and dimer models: 

(2) Zi(G):= Y, *(!*!)= E <M)=:Z£(r G ). 

eeZ 1 (G;Z 2 ) M6M(r G ) 
[?] = <* [M + M ]=a 

Remark. The construction described here is not the one given by Fisher 
in [7J. In his original construction, each vertex of degree n gives rise to 
3n — 6 vertices. It is therefore simpler for vertices of degree n < 5, in 
particular for square lattices. However, the construction presented here is 
simpler for vertices of degree n > 7, and overall more suited for the purpose 
of this paper. 

4.2. The geometric treatment of the dimer model. In this subsection, 
we recall the main features of the geometric method initiated in [3j , following 
the simplified version given in [2]. We refer to these articles for further 
details, and to [U [22| [9] for earlier treatments of the same problem. 

4.2.1. Dimers and Pfaffians: Kasteleyn theory. Let V be a finite graph. The 
aim is to compute the dimer partition function 

Z D (T) = x(M), where x(M) = ]J x e - 

MeM(r) eeM 

Kasteleyn's method is based on the following beautifully simple computa- 
tion. If there exists a perfect matching, then the number of vertices of T is 
even. Enumerate them by 1, 2, ... , 2n, and fix an arbitrary orientation K of 
the edges of T. 

Definition 2. The associated Kasteleyn matrix A K {T) = (Afj) is the 2n x 
2n skew-symmetric matrix whose coefficients are given by 

Afj = efj{e)x e , 

e 

where the sum is over all edges e in V between the vertices i and j, and 




1 if e is oriented by K from % to j; 
— 1 otherwise. 



Recall that the Pfaffian of a skew-symmetric matrix A = (A{j) of size 2n 
is given by 

Pf(A) = Y, (- 1 ) SgnW 4(l),(2) •••4,(2n-lM2n), 

Hen 
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where the sum is over the set II of perfect matchings of {1, ... , 2n} and a 
is a permutation of {1, . . . , 2n} representing the matching [a]. In the case 
of A (r), a matching of {1, ... , 2n} contributes to the Pfaffian if and only 
if it is realized by a perfect matching on T, and this contribution is ±x(M). 
More precisely, 

(3) Pi(A K (T))= £ e K (M)x(M), 

MeM(r) 

where the sign e K (M) can be computed as follows: if the perfect matching 
M is given by edges ei, . . . , e n matching vertices if and jg for I = 1, . . . , n, 
let a denote the permutation sending (1, . . . , 2n) to (h,ji, . . . ,i n ,jn), and 
set 

n 

(4) e K {M) = {-lf^J[sf dl {e i ). 

t=\ 

The problem of expressing Z D (G) as a Pfaffian now boils down to finding 
an orientation K of the edges of V such that e K (M) does not depend on M. 

Obviously, any perfect matching M can be considered as a cellular 1- 
chain M G Ci(T; Z2) such that <9M = t>, the sum being on all vertices of 
r. Hence, given any two perfect matchings M,M', their sum M + M' is a 
1-cycle. The connected components of this 1-cycle are disjoint simple loops 
of even length; let us denote them by {Cj}j. An easy computation shows 
that 

(5) £ K (M)e K (M') = H(-l) nK ^+\ 

i 

where n K (Ci) denotes the number of edges of C% where a fixed orientation 
of Ci differs from K. (Since Cj has even length, the parity of this number is 
independent of the orientation of Cj.) Therefore, we are now left with the 
problem of finding an orientation K of T such that, for any cycle C of even 
length such that T\C admits a perfect matching, n (C) is odd. Such an 
orientation is called a Pfaffian orientation. By the discussion above, if K is 
a Pfaffian orientation, then Z D (T) = \Pi{A K {T))\. 

Kasteleyn's early triumph was to prove that every planar graph admits a 
Pfaffian orientation. More precisely, let T be a graph embedded in the plane. 
Each face / of V C M? inherits the (say, counterclockwise) orientation of M 2 , 
so df can be oriented as the boundary of the oriented face /. 

Kasteleyn's Theorem ([14, 15J). Given Y C M 2 , there exists an orienta- 
tion KofT such that, for each face fofTc M 2 , n K (df ) is odd. Further- 
more, such an orientation is Pfaffian. 

An amazing consequence of this result is that it enables to compute the 
partition function of the dimer model on a planar graph in polynomial time. 
There is no hope to extend this result to the general case. Indeed, some 
graphs (such as the complete bipartite graph -^3,3) do not admit a Pfaffian 
orientation. More generally, enumerating the perfect matchings on a (bipar- 
tite) graph is a #P-complete problem [23]. It turns out that Kasteleyn's 
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method does extend to surfaces, but one needs to compute many Pfaffians. 
This is the aim of the following two paragraphs. 

4.2.2. Kasteleyn orientations as discrete spin structures. First note that any 
finite connected graph T can be embedded in a closed oriented connected 
surface £ as the 1-skeleton of a cellular decomposition X of S. (This simply 
means that the complement of T in X consists of open 2-discs.) In this 
paragraph, we shall use the same notation X for the graph embedded in the 
surface and for the induced cell complex. 

An orientation K of the 1-cells of such an embedded graph is called a 
Kasteleyn orientation on X if, for each 2-cell / of X, the following con- 
dition holds: the number n K (df) of edges in df where K disagrees with 
the orientation on df induced by the counterclockwise orientation on /, is 
odd. Given a Kasteleyn orientation on X, there is an obvious way to obtain 
another one: pick a vertex of X and flip the orientation of all the edges 
adjacent to it. Two Kasteleyn orientations are said to be equivalent if they 
can be related by such moves. Let us denote by %{X) the set of equivalence 
classes of Kasteleyn orientations on X. 

Proposition 4.2. A embedded graph X admits a Kasteleyn orientation if 
and only if X has an even number of vertices. In this case, the set %{X) is 
an affine H 1 ^;^) -space. 

The easy proof can be found in [31 Section 4]. 

Now, the game consists in trying to encode combinatorially a spin struc- 
ture on a surface S, or equivalently, a vector field on £ with isolated zeroes 
of even index. Let us begin by fixing a cellular decomposition J of S. 

• To construct a (unit length) vector field along the 0-skeleton X°, we 
just need to specify one tangent direction at each vertex of X. Such an 
information is given by a perfect matching M on X 1 : at each vertex, point 
in the direction of the adjacent dimer. 

• This vector field along X° extends to a unit vector field on X 1 , but 
not uniquely. Roughly speaking, it extends in two different natural ways 
along each edge of X 1 , depending on the sense of rotation of the resulting 
vector field. We shall encode this choice by an orientation K of the edges 
of X 1 , together with the following convention: moving along an oriented 
edge, the tangent vector first rotates counterclockwise until it points in the 
direction of the edge, then rotates clockwise until it points backwards, and 
finally rotates counterclockwise until it coincides with the tangent vector at 
the end vertex. This is illustrated in Figure 

• Each face of X being homeomorphic to a 2-disc, the unit vector field 
defined along X 1 naturally extends to a vector field on X, with one 
isolated zero in the interior of each face. One easily checks that for each 
face / of X, the index of the zero of Aj^ in / has the parity of n K (df) + 1 
(see [21 Lemma 4.1] for the proof). Therefore, this vector field defines a spin 
structure if and only if K is a Kasteleyn orientation. 

To summarize, a perfect matching M on X 1 and a Kasteleyn orientation 
K on X determine a spin structure on S. By Johnson's theorem (recall 
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Figure 3. Construction of the vector field along the 1- 
skeleton of X. 

Subsection 13. ip . it automatically defines a quadratic form : Z2) — > 

%2- This form is characterized by the following property: if C is an oriented 
simple closed curve on X 1 , then 

(6) qU[C])=n K (C)+t M {C) + l (mod 2), 

where £m(C) denotes the number of vertices in C such that the adjacent 
dimer of M points out to the left of C. (See [21 Lemma 3.3].) This fact 
easily implies the following correspondence theorem. 

Theorem 4.3. Let X be a cellular decomposition of an oriented closed sur- 
face S. Then, any perfect matching M G M(X X ) induces an H 1 ^;^)- 
equivariant bijection 

^ M ■ %(X) -> Q(S) = S(S), [K] 1 y q^ 

from the set of equivalence classes of Kasteleyn orientations on X onto the 
set of spin structures on S. □ 

4.2.3. The Pfaffian formula. As a direct consequence of Equation (|6j), we 
obtain "for free" the following non-trivial combinatorial result. 

Proposition 4.4. Let K be a Kasteleyn orientation on X , and M be a 

perfect matching on X 1 . Given a homology class a G i2i(£;Z2), represent 
it by oriented simple closed curves C\, . . . , C m in X 1 . Then, the equality 

m 

qM(u)=Y,( nK ( C i) +i M(Ci) + l)+ C i' C i ( mod2 ) 

1=1 l<i<j<m 

determines a well-defined quadratic form q^: H\(L;1,2) ~^ ^2- D 

We shall now use this combinatorial information, together with the results 
and notation of Paragraph 14.2. 1| to derive our Pfafhan formula. 

Let r be a finite connected graph. If T does not admit any perfect match- 
ing, then the partition function Z D (T) is obviously zero. So, let us assume 
that r admits a perfect matching Mq. Enumerate the vertices of V by 
1, 2, . . . , 2n and embed T in a closed orientable surface S of genus g as the 
1-skeleton of a cellular decomposition X of S. 

Since T has an even number of vertices, the set %(X) is an affine -ff 1 (S; Z2)- 
space (Proposition I4.2| ). For any Kasteleyn orientation K, the Pfaffian of 
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the associated weighted skew-adjacency matrix satisfies 

e K (M )Pf(A K (T)) i Y, e K (M )e K (M)x(M) 
MeM(r) 

§ £ (-l)^(«^)+i) x(M)) 
AfeM(r) 

where the Cj's are the connected components of the cycle M + Mq 6 
Ci(X;Z2). Note that given any vertex of Cj, the adjacent dimer of Mo lies 
on Ci, so that £M (Ci) = 0. Since the cycles C{ are disjoint, Proposition 14.41 
gives 

J> X (Q) + 1) = J>*(Ci) + e Mo (Ci) + 1) = <4([M + M ]). 

i i 

Therefore, every element [K] of %{X) induces a linear equation 

(7) e K (M )Pf(^(r)) = (-i)^° (a) ^a (r), 

ae/fi(E;Z 2 ) 

where Z^(T) = Y1[M+M ]=a X (M), the sum being over all M £ M(r) such 
that [M + Mo] = a. Equation (JTJ) now immediately yields the following 
Pfaffian formula. 

Theorem 4.5. Let T be a graph embedded in a closed oriented surface £ of 
genus g such that E\T consists of open 2-discs. Then, the partition function 
of the dimer model on V is given by the formula 

Z ^( r ) = i E (-i) Ari(q ^ ) e K (M )Pf(A K (T)), 
[K]eX(X) 

where the sum is taken over all equivalence classes of Kasteleyn orientations, 
and Arf(g) £ Z2 denotes the Arf invariant of the quadratic form q. □ 

4.3. Identification of the Kac-Ward and Fisher-Kasteleyn matrices. 

We are now ready to put together the results of the two previous subsections 
to obtain the geometric proof of Theorem 12.11 

Let G be a finite graph. First, number the vertices of G. Then, embed 
G in a closed oriented surface E, and blow up each vertex v G V{G) into a 
cluster as described in Figure[TJ This gives another graph Tq, also embedded 
in S. Enumerate the vertices of in each cluster by {a±, 61, ... , a n , b n }, as 
illustrated in Figure HI This, together with the numbering of the vertices of 
G, defines an ordering of the vertices of Tq. 

Let K be any Kasteleyn orientation on Tq C S. Replacing K by an 
equivalent Kasteleyn orientation, it may be assumed that it orients the edges 
in each cluster as described in Figure HI Let A k (Tg) denote the associated 
Kasteleyn matrix (recall Definition [2|) . 

Also, let \f, Io denote the vector field along T^r induced by K and by the 
preferred perfect matching Mo (Figure [2]) , as described in Paragraph 14.2.21 
This vector field extends to a nowhere vanishing vector field in the neighbor- 
hood of each cluster, thus defining a vector field A along G. Let I — T X (G) 
be the associated generalized Kac-Ward matrix (recall Definition [1]) . 
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Figure 4. The numbering of the vertices, and a Kasteleyn 
orientation at a cluster. 

Proposition 4.6. We have the equality Pf(A K (T G )) = det(J - T X (G)) 1 / 2 . 

Note that with our enumeration convention, the sign e k (Mq) is equal to 
+1 (recall Equation (jl|)). Therefore, we immediately obtain Theorem l2.lt 

M}_ £ (-i) Arf (^o)pf(^ ( r G )) 

[K]£X(X) 

m h E (-i) Arf(A) Pf(^(r G )) 

Aes(s) 

IPJL £ (-l) Arf Wdet(J-T A (G)) 1/2 . 
Aes(s) 

Proof of Proposition \4-6\ First note that the constant coefficient of the poly- 
nomial det(I — T A (G)) 1 / 2 is (by choice of the square root) equal to +1. 
On the other hand, the constant coefficient of Pf(A K (Tg)) is equal to 
rLeV(G) (Ateg(«))i where A n is the 2n x 2n Kasteleyn matrix associated to 
a cluster coming from a vertex f £ V(G) of degree n. One checks by direct 
computation that Pi(A n ) is equal to +1 for any n. Hence, both polynomials 
in the statement of Proposition 14.61 have the same constant coefficient. It is 
therefore enough to show that their squares coincide, that is, to check the 
equality 

(*) det(^(r G )) = det(/-T A (G)). 

Then, observe that both sides of this equation are left unchanged when 
adding an edge e to G C S with associated weight x e = 0. Indeed, if G' C S 
denotes this new graph, then I—T X (G') = (I —T x (G))®l2, where I2 denotes 
the size 2 identity matrix. As the left-hand side of (*) is equal to 

det(^(r G )) i (£(-l)<> ( ^f (F G )) 2 i (£(-l)^o^(G)) 2 , 

a a 

it is obviously also left unchanged. Therefore, it may be assumed without 
loss of generality that all vertices v S V(G) have even degree. (This as- 
sumption is by no mean crucial, but it will simplify the computations to 
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come.) Indeed, given two vertices v,v' G V(G) of odd degree, double each 
edge in a path from v to v'. The new graph will now have even degree at 
both v and v', while all the other degree parities are unchanged. As the sum 
J2veV(G) deg(v) = 2|J5(G)| is even, the number of vertices in G with odd 
degree is even, and the procedure above allows to get rid of all of them. 

Next, it may be assumed that G has no loop. (Again, this is only for 
notational simplicity in the computations to come.) Indeed, given a loop e 
at a vertex v S V(G), add a vertex v' in it: this divides e into two edges 
to which we assign the weights x e and 1. Let G" C £ denote this new 
graph. It is an amusing exercise to check that the right-hand side of (*) is 
unchanged when replacing G by G": details are left to the reader. As Z l a (G) 
and Z^(G") obviously coincide, the invariance of the left-hand side of (*) 
follows directly from the equation displayed above. 

So, let G be a finite graph with no loops, and with each vertex of even 
degree. The strategy to prove (*) will be to transform A k (Tq) into the 
matrix I -T X {G) using moves that do not change the determinant. Let us 
focus on an arbitrary vertex v £ V(G) of degree n. By definition of K and 
numbering convention of the vertices in the corresponding cluster (recall 
Figure HD, and since there is no loop at v, the corresponding part of the 
matrix A k (Tq) has the following form 

\ 
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where rj, r. 



(resp. 



denote rows (resp. columns) of the matrix. Using 
elementary row and column operations and the fact that n is even, one can 
easily modify A (Tq) to obtain 
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It is of historical interest to note that this is precisely the type of matrices 
first considered by Kasteleyn [13] and Fisher [6] to solve the Ising model 
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on square lattices (case n = 4). It is not strictly speaking a Kasteleyn 
matrix, as the corresponding graph T' G does not embed in the same surface 
as G anymore: every vertex v G V(G) of degree n gives rise to a complete 
graph on n vertices. However, its Pfaffian still gives the right answer for the 
Ising model on G. It is only later that Fisher [2j found his correspondence 
G i — y Tq where Tq remains embedded in the same surface as G. 

Let L n denote the n x n skew-symmetric matrix appearing above. Its 
determinant is equal to 1, and its inverse is given by 





(0 


-1 


-1 . . 


• -!\ 




1 





-1 . . 


. -1 




1 


1 


. . 


. -1 




V 1 


1 


1 . . 


0/ 



Therefore, 

det(^(r G )) = det(A K (T' G )) = det <A K (V G ) ■ L^). 

veV(G) 

Let B K denote this latter matrix. Its rows and columns are indexed by 
V(T' G ) = Uuey(G) C(v) C V(Tg), where C(v) denotes the set of vertices 
{ai,a2 . . . ,adeg(i>)} i n t ne cluster associated to v G V(G) (recall Figured]). 
The coefficient of B K corresponding to a G C(v) and a' G C(v') is equal to 



B 



1 if a = a' , 

£ aa"( e )^a"a' x e if there is an a" G C(y') joined to a in Tq, 
otherwise, 

where e G E(Fq) denotes the edge joining a and a", e^«(e) = 1 if e is 
oriented by K from a to a", e^/(e) = —1 otherwise, and 

+1 if a' < a" in the linear ordering of the vertices of C(v'), 
-1 if a' > a", 
if a' = a". 

Finally, let B K denote the matrix obtained from B K by multiplying by 
i = (resp. by —i) each row (resp. each column) corresponding to a 

vertex a G V(Tq) whose adjacent edge in E(G) C E(Tg) is oriented by K 
towards a. We shall write K{a) = 1 if this condition holds, and K(a) = 



otherwise, so that e 



^„(e) = (-l) K ( a ). Obviously, 
det(B K ) = det(B K ) = det(A K (T G )), 
and the coefficient of B K corresponding to a G C{v) and a 1 G C(v') is 
1 if a = a' , 

(-j)*(°)+*( a ')4,„ a , x e if there is an a" G C(t/) joined to a in r G , 
otherwise. 



B. 



K 



Now, recall that the matrix / — T A (G) is indexed by the set E of oriented 
edges of G. Let (p: E — > V(T' G ) denote the bijection mapping each oriented 
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edge e to its starting vertex in Fq- By definition of the matrix T X (G), and 
since G has no loops, 



(7 - T\G)) t 



1 if e = e', 

— exp (|a^(e, e')) x e if /(e) = s(e') but e ; ^ — e, 

otherwise. 



For the two oriented edges e, e' of G, let s(e) = v, s(e') = u ; denote their 
starting vertices in G, and (p{e) = a, y(e') = a' their starting vertices in Tq- 
Then, the condition /(e) = s(e') and e' ^ — e is equivalent to the condition: 
the oriented edge e of Tq joins a to some a" E C(t/) with a" 7^ a'. Therefore, 
the coefficient (I — T x (G)) e ^ e i is non-zero if and only if B^,-, ,~ is non-zero. 

Finally, using the construction of the vector field A = A^ o induced by K 
and Mo, one easily checks that for any e, e' E E satisfying this condition, 

-exp(§a A (e,e')) = (-i)^)+^'H a „ a ,. 

Therefore, the matrices I — T X (G) and B K coincide. This concludes the 
proof of Proposition 14.61 □ 

As mentioned in the remark at the end of Subsection 14. 1|. the correspon- 
dence G h- > Tq used here is a variation of Fisher's original construction [TJ, 
that could have been used just as well. Assuming that all the vertices of G 
are of even degree, there is yet another variation of the Fisher correspon- 
dence where each vertex of degree n gives rise to 2n — 2 vertices (see [1]). 
It is also possible to use this version of the correspondence: the proof of 
Proposition 14.61 only requires minor adjustments to be adapted to such vari- 
ations on the theme. This yields an identification of the Kac-Ward method 
with the particular version of the Pfaffian method described in [3]. Such a 
result is by no mean obvious: the special case of the planar square lattice 
alone is quite non-trivial, being the subject of the article [25] , 
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